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Abstract. We have studied the Dalitz plot of the e+e" — » -k^tt^j events collected at -^/s ~ with the 
KLOE detector. In the dipion invariant mass (M,r7r) region below 700 MeV, the process under study is 
dominated by the non-resonant process e'^e~ — > wtt" with oj — > n'^y whereas, for higher Mt^t, values, the 
radiative <f) decay to the /o(980) is the dominant mechanism. Different theoretical models are used to fit 
the Dalitz plot, taking also into account a possible contribution of the cr(600). For each model, we extract 
the /o(980) mass and its coupling to tttt, KK and to the 4>. 
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1 Introduction 

Interest in light scalar mesons remains intense in hadron 
spectroscopy due to a lack of elucidation on their nature. 
There is a possibility that some of them arc, in fact, ex- 
otic particles. There are several models to describe their 
structure, such as ordinary qq mesons, qqqq states or KK 
molecules [1-3]. Operating at the e+e~ Prascati (^factory 
DA<?NE [4], the KLOE experiment [5] is ideally suited 
for the study of these particles, since the radiative de- 
cays of the ^ into two pseudoscalar mesons is dominated 
by a scalar meson (S) exchange in the intermediate state 
(0 5*7 TTTTj/riTTj/KKj). For the 7r°7r°7 final state, 
the possible scalar contributions are from the well estab- 
lished /o(980) and from the more controversial (t(600), 
purportedly observed by the E791 and BES collabora- 
tions [6,7]. The non-resonant e+e^ wtt*^ 7r°7r"7 re- 
action also contributes to the same final state. 



The 7r°7r°7 final state had been studied at KLOE using 
16 pb~^ of 2000 data [8]. The resulting ratio between the 
foKK and /otttt couplings, together with the large value 
of the BR, favoured the qqqq composition of the /o(980), 
while the shape of the tt^tt^ invariant mass suggested a 
possible contribution also from cr(600). In the present pa- 
per, the analysis is repeated with a statistics about thirty 
times larger, thus allowing us to study this reaction in 
much greater detail. A common set of cuts and algorithms 
for the resonant and non-resonant processes has been de- 
veloped so that, by fitting the Dalitz plot, the differential 
cross sections of the two components are extracted. A de- 
tailed technical description of this analysis is in Refs. [9, 
10]. 

2 Experimental setup 

Data were collected with the KLOE detector at DA^PNE, 
the Frascati e+e~ ^factory, which operates at a center 



of mass energy y/s = ^ 1020 MeV. The beams col- 
lide with a crossing angle of (tt — 25) mrad, producing 
<p mesons with a smah momentum {p^ ~ 13 MeV/c) in 
the horizontal plane. The KLOE detector (see Fig. 1) is 
inserted in a 0.52 T magnetic field. It consists of a 2 m 
radius drift chamber (DC) [11], with full stereo geometry 
using helium based gas mixture, surrounded by a fine sam- 
pling lead/scintillating fibers electromagnetic calorimeter 
(EMC) [12], divided into a barrel and two endcaps, with 
a hermetic coverage (98% of the solid angle) and a very 
high efficiency for low energy photons. Since the channel 
e^e^ — > 7r'^7r''7 under study is fully neutral, its analysis is 
based mainly on the EMC performance. The arrival times 
of particles and the positions in three dimensions of the 
energy deposits are obtained from the signals collected at 
the two ends of the calorimeter modules, with a granular- 
ity of ~ (4.4 X 4.4) cm^, for a total of 2440 cells arranged 
in five layers. Cells close in time and space are grouped 
into a calorimeter cluster. The probability of a photon to 
fragment in more than a cluster (splitting) is reduced by 
employing a special recovery algorithm. The cluster energy 
E is the sum of the cell energies, while the cluster time T 
and its position R are energy weighted averages. Photon 



energy and time resolutions are ue/E — 5.7%/ ^E (GeV) 



and ax = 57 ps/^£: (GeV) e 100 ps, respectively. The 
KLOE trigger [13] is based on the detection of two energy 
deposits (called sectors) with E > MeV for barrel and 
ii^ > 150 MeV for endcaps. Events with only two fired 
trigger sectors in the same endcap are rejected, being this 
topology dominated by machine background. Recognition 




-6 m- 



Fig. 1. Vertical cross section of the KLOE detector. 

and rejection of cosmic-ray events is also performed at 
the trigger level, selecting events with two energy deposits 
above a 30 MeV threshold in the outermost calorimeter 
layer. Moreover, to reject residual cosmic rays and ma- 
chine background events, an offline software filter uses 
calorimeter and DC information before track reconstruc- 
tion [14]. 

The machine parameters (center of mass energy -y/s, 
4> momentum and beams interaction point) are measured 
on-line from the analysis of Bhabha scattering events in 
the barrel. The average value of the center-of-mass energy 
is evaluated with a precision of 30 keV per run, each cor- 
responding to ~ 100 nb^^ of integrated luminosity. To 
calibrate the absolute beam energy scale we fit the line 
shape of ^ 777 events [9] . Comparing the value obtained 
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for with the precise measurement of the CMD-2 ex- 
periment obtained with the depolarization method [15], a 
shift of +150 keV is found and is corrected in our analysis 
accordingly. 

Prompt photons are identified as neutral particles with 
/3 = 1 originated at the interaction point, by requiring 
\T — R/c\ < min(5crT,2 ns), where T is the photon flight 
time and R the corresponding path length. The photon 
detection efficiency is 90% for E^=20 MeV, and reaches 
100% above 70 MeV. The sample selected by the timing 
requirement has loss than 0.6% contamination per event 
due to accidental clusters from machine background. 

3 Event selection 

All the available statistics collected in 2001-2002 data- 
taking periods, corresponding to 450 pb~^, has been an- 
alyzed by grouping all runs in center of mass energy bins 
of 100 keV. This was done to take into account the ~ 2 
MeV spread in the center of mass energy present in the 
data set. For this analysis, only those runs belonging to 
the bin with the highest statistics have been used for fit- 
ting the Dalitz plot. This sample corresponds to 145 pb~^ 
collected at ^ = (1019.7+1019.8) MeV. 

The response of the detector to the decays of interest 
was studied by using the KLOE Monte Carlo (MC) sim- 
ulation program [14] . The MC takes into account changes 
in the machine operation and background conditions, in 
order to reproduce real data on a run-by-run basis. For 
the present analysis, an MC sample for both signal and 
backgrounds is produced. The corresponding integrated 



luminosity is five times that of the collected data, ex- 
cept for e+e~ — > 77 events that are produced at a 1:1 
rate. For the simulation of signal events, the M^^r spec- 
trum for the (j) ^ Sj ^ wwy (S-y) process is produced 
according to the shape obtained from 2000 data, while 
the e+e~ — > W7r° — > tt^tt^j (wtt) generator is based on the 
Vector Meson Dominance (VMD) description of the three 
body decay according to Ref. [16]. 

The data analysis consists of four steps: 

1. an acceptance selection of five prompt photons with 

>7 MeV and a polar angle satisfying the require- 
ment I cos9^\ < 0.92; 

2. a kinematic fit (Fitl) imposing total 4- momentum con- 
servation; 

3. a pairing procedure of photons to 7r°'s, where the pho- 
ton combination minimizing a pseudo-x^ built using 
the invariant mass of the two 77 pairs, Xpair> se- 
lected as the good one; 

4. a second kinematic fit (Fit2), where the constraints on 
the 7r° masses are also imposed. The selected events 
must then satisfy the requirements XFit2/^'^^ — ^ ^^'^ 
AMyj = \Myj — Mjrl < 5 (T.^,,, where M^y and a^-y are 
evaluated using the photon momenta from Fitl. 

A further cut is applied to reject the background from 
the e+e~ 77 process, which has a cross section much 
larger than the signal and where the three additional prompt 
photons could be generated either by radiation, by cluster 
splitting or by accidental coincidence with machine back- 
ground clusters. Such process is hugely reduced without 
losing efficiency for the signal, by rejecting events where 
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Fig. 2. Data-MC comparison after acceptance selection: nor- 
malized 'x from the first kinematic fit (top) and the minimum 
value of the pseudo-x^ used to pair photons (bottom). 

the energy sum of the two most energetic clusters in the 
event is greater than 900 MeV. 

The overall analysis efficiency for the identification of 
the signal is evaluated by applying the whole analysis 
chain to the and cjtt MC events: £57 — (50.3 ± 0.1)%, 
SujTx — (53.12 ± 0.05)%. The small difference is due to the 
characteristic energy and angular distributions of photons 
in the two kinds of events. After acceptance selection we 
start with a sample of 243,904 events of which 86,449 sur- 
vive the complete analysis chain. As shown in Figs. 2, 3 
excellent data-MC agreement is found, both after accep- 
tance selection and after applying the complete analysis 
chain. 

The background channels for the -k^-k^^ final state are 
listed in Tab. 1 together with the analysis efficiency and 
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Fig. 3. Top: normalized x'^ of the second kinematic fit after 
acceptance cuts. Bottom: normalized AM-,j of the selected 77 
pairs after XFit2 cut in linear (left) and logarithmic (right) 
scale. 

the corresponding signal to background ratio, before and 
after the application of the analysis cuts, evaluated using 
branching ratios (BRs) from Ref. [17] or KLOE measure- 
ments whenever available [18,19]. The (f> 7777^7 channel 
(rinj) is the only one that has the same five photon final 
state as the signal. The 777 —> 7r''7r''7r'^7 (7777) mim- 
ics the signal when there are lost or merged photons. The 
three photon final states {(f) — s- r/7 [7773], (j) — s- 7r'^7 [7r7] and 
e+e^ 77(7) [77]) produce five clusters due to split- 
ting or accidental coincidence with clusters produced by 
machine background. We have used background-enriched 
distributions to check the absolute yields and the Monte 
Carlo shapes. Each data distribution (Hdata) has been fit 
with two MC components: the background under consid- 
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eration (Hbckg) and all the other contributions, including 
the 7r°7r°7 signal (Hothcrs): Hdata = "i Hbckg + "2 Hothcrs- 
These distributions, used to fit the main background com- 
ponents after applying the ai scale factors, are shown in 
Fig. 4. The dominant contribution 7777 is verified by study- 
ing the events in the region 4 < Xp^^^/'^'^^ < 20. For all 
other background sources, we keep the standard analysis 
cuts and build a specific i minimizing the difference be- 
tween the reconstructed and the true mass of the interme- 
diate particles (77 and tt^) in the corresponding hypothesis. 
The values of ai obtained for all the backgrounds are listed 
in Tab. 1 together with the uncertainties on the BRs. We 
do not apply these scale factors, but we use them in the 
evaluation of the fit systematics as discussed in Par. 4.3. 
Note that for the dominant contribution (7777) a scale fac- 
tor statistically consistent with the expected rate is found. 

In order to check the relative contribution of the two 
signal processes, their angular distributions are studied. 
Both the photon polar angle (9) and the minimum angle 
between the photon and the 7r°'s in the tt^tt" rest frame 
(ip) are expected to show a different behaviour due to 
the spin of the intermediate particles involved. To first 
order the interference between the two processes can be 
neglected and they can be separated on event by event ba- 
sis by looking at the mass of the intermediate state. The 
Lo mass is reconstructed by selecting the best match of 
the two TT^j combinations. After background subtraction, 
events with \M^^ — M^^l < 3aM^^ are classified as cjtt, 

^ <^Mu, = 9-5 MeV is the convolution between the experimen- 
tal resolution and the u width. 




5 10 15 20 5 10 15 20 

2 2 

X % 

Fig. 4. Distributions used to check background contami- 
nation. Top-left: photon energy distribution for events with 
4 < Xpit2/iidf < 20. Minimum value of the used to pair pho- 
tons in the 7773 (top-right), 717 (bottom- left) and rj-K'y (bottom- 
right) hypotheses. The 87 distributions include 7/73, 717 and 77 
processes. 

while all the rest is called 6*7. In Fig. 5 the data-MC com- 
parison for the cos 9 and cos ip angular distributions are 
shown for both 6*7 and lutt processes. The simple super- 
position of the 5*7 and lutt MC shapes fits rather well the 
data, suggesting that the contribution of the interference 
term is small. 

The two kinematic variables chosen for the construc- 
tion of the Dalitz plot are the invariant masses of the 
two tt'^'s, A/tttt, and of the two possible 7r°7 combinations, 
M^-y. We have therefore two entries per event. The binning 
choice, driven by the mass resolutions obtained by Monte 
Carlo for the signal, is 10 MeV for and 12.5 MeV for 
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Table 1. Analysis efficiency for background events. The signal over background ratio before (S/B) and after (S/Bana) the 
application of the analysis cuts is also shown. In the last two columns the error on the expected rate due to the uncertainties 
on the used BRs and the scale factor obtained from the fit to the background enriched distributions are reported. 



Process 


£ ana 




S/B 


S/Bana 


(5BR/BR 


Ql 


(j> — > rj-K'^'y 777r"7 


(23.2 ±0.1 


)% 


8.5 


19 


9.5% 


0.86 ±0.02 


(j> ^ rij 7r*V'V"7 


(8.51 ±0.02) ■ 


10"^ 


0.06 


4 


3.0% 


1.064 ± 0.002 


(p ~* Vl ^ 777 


(8.15 ±0.05) • 


10-'' 


0.06 


30 


3.5% 


0.86 ±0.02 


(j> n'^'y 


(3.07 ±0.06) ■ 


10-'' 


0.2 


350 


10.0% 


2.35 ±0.02 


e+e~ 77 


(0.19 ±0.01) • 


10-^ 


0.002 


400 




1.85 ±0.03 
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Fig. 5. Data-MC comparison for con (top) and 5*7 (bottom) 
events selected by looking at the mass of the intermediate state 
and assuming the interference term to be negligible. Left pan- 
els: cos distributions. Right panels: cos for the primary pho- 
ton. 

Afjr-y. The data density is shown in Fig. 6 before and after 
background subtraction. The two projections arc shown in 
Fig. 7. After background subtraction the number of events 
in the Dalitz plot is 128, 529 ± 659. 



Fig. 6. Dalitz plot in logarithmic scale before (left) and after 
(right) background subtraction. The two bands in the region 



M^^ < 700 MeV are due to e+e" 



LOTv" events. 



The analysis efficiency as a function of Af^jr and M^^^ 
is evaluated by Monte Carlo, with corrections based on 
data control samples. In Fig. 8 the dependence of the se- 
lection efficiency on M^^t^ and Mtt-/ is shown for the 57 and 
WTT final states. Both processes exhibit a rather flat depen- 
dence on the Dalitz plot variables. The different shape is 
related to their angular distribution and to their different 
probability of producing photons from initial state radi- 
ation (ISR). The main source of data-MC differences is 
due to the photon detection efficiency, which is measured 
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Fig. 7. Dalitz plot projections in M^^ (top) and M^-, (bottom) 
variables before background subtraction. 

and applied, as a correction, to the Monte Carlo. The MC 
trigger, cosmic ray veto and event classification filter ef- 
ficiencies are checked using prescalcd data samples. The 
overall correction factor, i?soi = f .022 ± 0.004, is applied 
to the MC analysis efficiency. 

4 Fit to the Dalitz plot 
4.1 Fitting method 

The double differential 7r°7r°7 cross section is written as 
the sum of the scalar term, the VMD contribution and of 
their relative interference as follows: 



da 



(1) 



da"- 



da^ 



dM^T^dAIj^^ 



± 



da' 



dMTTTrdM^j 



Fig. 8. Analysis efficiency for 5*7 (empty triangles) and ujtt 
(black dots) MC events as a function of AIt„t (top) and M,r7 
(bottom). 

The suffixes S and V stand for 57 and VMD terms while 
/ represents the interference. To fit the Dalitz plot data 
density, the theory is folded with the reconstruction efh- 
ciencies of the two processes and with the probability for 
an event to migrate from a Dalitz plot bin to another one. 
The expected number of events for a given reconstructed 
Af^jr, Mttj bin i, N^^^, is then computed from the total 
integrated luminosity, -Lint, as follows: 

Nr = imt X (2) 

X! (/j^ + -^Z ^5 ^3 + ■^Z • 

where fj is the integration of the differential cross section 
evaluated in the j*'' bin, including the effect of the ISR. 
For each bin, Sj is the analysis efficiency while Aij is the 
smearing matrix, representing the probability for the sig- 
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pair 

Fig. 9. Top: Monte Carlo Zlxpair shapes for events with correct 
and wrong photon pairing. These two templates are used to 
fit the z^Xpair data distribution (bottom). A better data-MC 
agreement is obtained after the fit. 



nal event to migrate from the j to the i bin, either due 
to resolution or to a wrong reconstruction of the event. 
This matrix has been evaluated by Monte Carlo: about 
85% of the events are on the diagonal or close to it, within 
±1 bin. The remaining 15% is mostly due to events where 
photons are incorrectly paired to tt^'s. A dedicated data- 
MC comparison has been performed to calibrate the frac- 
tion of good/bad pairing. The difference between the two 
lowest values of pseudo-x^ used to pair photons, zixpair' 
is fit with a superposition of the two templates obtained 
for MC events with good and bad photon pairing (Fig. 9). 
A data-MC difference of 1.08 ± 0.02 is found and is taken 
into account in evaluating the systematic error. 



4.2 Theoretical models 

A more explicit formulation of the differential cross section 
(1) is reported in Appendix A. Concerning the scalar term, 
we use two different approaches for the description of the 
amplitude: the Kaon Loop model (KL) [20-22], where the 
scalar is coupled to the (j) through a charged kaon loop, and 
another formulation called No Structure (NS) [23], where 
the (jjSj coupling is point-like and the scalar is described 
as a Breit-Wigner with a mass dependent width added to 
a complex function, to allow an appropriate behaviour of 
the resulting shape at low masses. 

4.2.1 Kaon Loop model 

As scalar term of the Kaon Loop model we use the one 
described in Ref. [22], where the amplitude of the signal 

M57 = e'^^g{M^^) i^ gsK+K- Ggl, gs'^o.^ , (3) 

implies the mixing of two scalar states, namely the /o(980) 
and the cr(600), represented by the Gss' matrix. The g{M-^^) 
function describes the kaon loop and gsTnr/KK a-re the cou- 
plings of the scalars to the rnr/KK mesons. Differently 
from the past, where the phase Sb took into account only 
the elastic scattering due to tttt, in this new formulation 
the scattering due to KK is also considered. 

Using this theoretical framework, our tt^tt^ mass spec- 
trum obtained from 2000 data [8] , has been fit by the au- 
thors of Ref. [22] together with tttt scattering data [24-28] , 
providing ten sets of parameters which are able to describe 
both data samples. For all of them, the (7(600) coupling to 
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ii'ii' is small with respect to the one of TTTT and its mass lies fl'(/i/o7, and a subsequent coupling of the /o(980) to the 
in the 500-700 MeV range. This means a very broad me- tttt pair, g/gTrn- In the same form factor the possible non- 
son width, between 240 and 490 MeV. All the odd (even) resonant continuum background is also added as a series 
variants have a positive (negative) Qan+n- /dfoK+K- ratio, expansion in M„„ as follows: 

When fitting the Dalitz plot of our new data set, we ^scai lun" j^'' IJojui 

can not leave all /o(980) and (j(600) parameters free, as ^'^ fo~^ ' ^'^ 

the description of the clastic background and the cou- '^'^ '° ' 

plings and masses of the cr(600) meson are closely re- where, in the most general case, the background param- 

lated. We therefore proceed by using as free parameters eters are complex numbers and w^(M„^) is proportional to 

only the VMD ones and the mass and the coupling to pion momentum in the scalar rest frame, f ^ = ^Ml^jA - M^. 

7r+7r~ and KK of the /o(980), using the isospin rela- The propagator for the /o(980) resonance is described by 

tions gfoTr°w° = l/v^^/oTT+Ti— = l/v^ff/oTTTT aud QfoKK = a simple Breit-Wigner shape corrected by the Flatte con- 

9foK+K- = 9fgK0K°- ^'^^ cr(600) and the elastic tttt dition on the tttt and KK thresholds, i.e.: 

and KK scattering we use the values of Ref. [22] by re- y (j^ \ 

rf,{M^^) = g%,, + (5) 

peating the fit for all the ten available sets of parameters. '^'^ 
We obtain P(x^) ranging from 10"^ to 0.145; only the six ^SoKK SttM^^^ 

results with P(x^) > 1% are kept for the present discus- where v„^k,ko are complex numbers with an analytical 

sion. continuation under threshold and the coupling to tttt and 

Furthermore, the case of a single scalar contribution, KK have the same meaning as in the KL description. In 

the /o(980), in the Kaon Loop description has also been this model the fit parameters are Mf^, 5/0^+^-, gfoKK, 

tried by using the old KL parametrization [21]. The result- g,pfa-y and the background parameters ao, ai and bi. The 

ing P(x^) of the fit is not acceptable, showing the need of bo phase is fully determined as a function of the other 

introducing the cr(600). parameters. 



4.2.2 No Structure model 

In the NS description, the amplitude of the scalar term is 

proportional to a scalar form factor, Fq'^'^^, which is derived 

by assuming a direct coupling of the to the /o(980), 
^ The accepted results are the variants 1, 2, 3, 8, 9, 10 of 

Ref. [22]. 



4.3 Fit Systematics 

There are different sources of systematics affecting this 
analysis which can give rise to variations of the fit results. 
We describe here the most important ones. For each of 

them, the fit is repeated after changing the relevant quan- 
tity within its range of uncertainty. In Tabs. 2, 3 we show 
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Table 2. Fractional systematic error on fit parameters for the KL scalar term. Only variations above 0.1% are reported. 



Parameter 


Mfo 


9foK+K- 


5/o^+7r- 


Source 


Fractional systematic error 


Normalization 




±1.6% 


±0.7% 


Beam energy 




+2.4% 


-1.4% 


Photon efficiency 




-2.1% 


+ 1.4% 


cut 




-0.8% 




Smearing matrix 


-0.1% 




+1.4% 


Interference 




-1.1% 


+0.7% 


Background 


+0.1% 


+4.0% 


-3.5% 



Table 3. Fractional systematic error on fit parameters for the NS scalar term. Only variations above 0.1% are reported. 



Parameter 


Ml, 


9foK+K- 


f/oT+TT- 




ao 


ai 


6i 


Source 


Fractional systematic error 


Normalization 


±0.1% 


+ 155% 
- 50% 


+6.9% 

-1.5% 


+11.9% 

- 1.9% 


+34.4% 
- 5.5% 


+122% 
- 16% 


+57.3% 
-16.7% 


Beam energy 


-0.2% 


±82.5% 


+0.8% 


±1.9% 






- 3.3% 


Photon efficiency 


±0.2% 


-72.5% 


+3.1% 


±3.4% 


+17.9% 


+49.6% 


+29.3% 


cut 


-0.2% 


±57.5% 


+2.3% 


-2.3% 


+ 1.7% 


+17.6% 


+20.0% 


Smearing matrix 


-0.2% 


- 7.5% 


-2.3% 


-3.1% 


-10.0% 


-26.0% 


-28.0% 


Interference 




±82.5% 


+3.1% 


+5.4% 


+18.1% 


+69.5% 


+46.0% 


Background 


-0.4% 


+92.5% 




-2.3% 


-8.4% 


-8.4% 


+ 4.0% 



the corresponding percentage variation of the interesting 
free parameters for the KL, NS model respectively. 

1. Normalization 

The first effect considered is the normalization scale of 
the fit estimate on the event counting. When evaluat- 
ing N^^^ , two experimentally determined constants are 
used: the integrated luminosity and the value of the (j) 
leptonic width. The luminosity is known with a total 
error of 0.6% [29] while the leptonic width has been 
measured by KLOE with a 1.7% uncertainty [30]. The 



fit has been repeated by changing the value of both 
quantities of ilcr. 
2. Beam energy scale 

The beam energy scale also affects the fit due to the ex- 
plicit ^ dependence of the theoretical function. Bhabha 
scattering events allow a relative calibration of the en- 
ergy scale with a precision of 30 keV each 100 nb~^ 
of integrated luminosity. As mentioned in the intro- 
duction, an absolute calibration of 150 keV has been 
applied to match the measured value of M^. The fit 
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to the Dalitz plot has been repeated without applying 
this correction. 

3. Photon efficiency 

The data-MC correction of the cluster efficiency curve 
modifies the shape of the Dalitz plot. A different pa- 
rametrization of the cluster efficiency curves in the 
Monte Carlo [31] has been used to evaluate anew the 
smearing matrix, the analysis efficiency and the back- 
ground contribution. The fit has been repeated in these 
conditions. 

4. cut 

To test the systematic contribution of the chosen 
cut, we have repeated the whole analysis hardening the 
XFit2/^dof cut from 5 to 3. In this way the event count- 
ing is improved due to the large reduction of back- 
ground while the analysis efficiency is not as flat as 
before along the Dalitz plot due to the ISR tails. 

5. Smearing matrix 

From the measured quantity of wrong photon pairing 
(Sec. 4.1), the fraction of off-diagonal events in the 
smearing matrix has been increased by 8%. 

6. Interference 

In the standard fit function, the radiative corrections, 
the analysis efficiency and the smearing matrix used in 
the interference term are obtained from an MC sample 
of UTT events. The fit is repeated by using the corre- 
sponding quantities estimated with ^7 events. 

7. Background 

The scale factors ai, obtained when fitting the back- 



ground-enriched distributions, are applied to the resid- 
ual background contamination. 

As shown in Tabs. 2, 3, the KL fit provides stable re- 
sults when the systematic changes are applied. On the 
other hand, the NS fit shows very large parameter vari- 
ations due to the presence of the continuum background 
term. 

4.4 Fit results 

4.4.1 Kaon Loop model 

The Kaon Loop fit results are listed in Tabs. 4, 5 for the 
six accepted variants. In Fig. 10 the distributions of the 
data points for all slices of the Dalitz plot with the super- 
imposed fit function of the variant with the best P(x^) 
is shown. Here the M^^ vs M„„ histogram is sliced in 
M„„ and for each slice the M,^^ projection of the allowed 
phase-space region is plotted, one after the other. To un- 
derstand the relative importance of the different fitting 
terms, their contributions are shown in Fig. 12. left. As 
expected, the VMD (57) term is dominating in the region 
below (above) 700 MeV. The interference term is concen- 
trated around 600 MeV. As a fit result, we use the central 
value and the errors from the best fit, adding the system- 
atic error discussed in the previous section and an extra 
error associated with the theoretical model. This last error 
is cvahiatcd as the maximum variation between the cen- 
tral value obtained by the best fit (-P(x^) = 0.145) and 
the other five accepted fits. The extracted parameters are: 

Mf, = (976.8 ± 0.3 fit t°o% syst + 10.1 n,od) MeV 
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Fig. 10. Fit result for Kaon Loop model. Black dots are data while the sohd line represents the KL resulting shape. 



Table 4. Fit results for the scalar term in the Kaon Loop model. Different rows are obtained for the accepted sets of parameters 
from Ref. [22]. 



Fit 


Mto (MeV) 


9foK+K~ (GeV) 


(GeV) 




xVNdof 


Fix') 


Kl 


976.8 ± 0.3 


3.76 ± 0.04 


-1.43 ±0.01 


462 , 286 


2754/2676 


0.145 


K2 


986.2 ±0.3 


3.87 ±0.08 


-2.03 ±0.02 


485 , 240 


2792/2676 


0.058 


K3 


985.2 ±0.2 


4.92 ±0.06 


-1.92 ±0.01 


472 , 320 


2809/2676 


0.036 


K4 


982.3 ±0.4 


4.02 ± 0.06 


-1.76 ±0.02 


415,260 


2787/2676 


0.066 


K5 


983.3 ±0.2 


3.75 ± 0.02 


-1.40 ±0.01 


529 , 366 


2823/2676 


0.024 


K6 


986.9 ±0.1 


3.28 ±0.05 


-1.90 ±0.01 


566 , 264 


2799/2676 


0.048 
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Table 5. Fit results for the VMD parametrization. The first six rows are obtained in the KL approach for the accepted sets of 
parameters from Ref. [22] while the last line is the result of the NS fit. 



Fit 


Ctpn 


C^^ (GeV"'') 








K n 


(McV) 


Kl 


0.58 ±0.11 


0.850 ±0.010 


0.46 ±0.13 


0.260 ±0.185 


3.11 ±3.12 


33.0 ±9.7 


782.52 ± 0.29 


K2 


0.68 ±0.03 


0.832 ±0.003 


0.30 ±0.05 


0.061 ±0.211 


3.14 ±3.08 


23.6 ±4.1 


782.20 ±0.11 


K3 


0.66 ±0.17 


0.836 ± 0.004 


0.33 ± 0.08 


0.084 ± 0.056 


3.14 ±3.14 


25.2 ± 6.2 


782.26 ± 0.28 


K4 


0.64 ± 0.05 


0.836 ± 0.002 


0.38 ± 0.06 


0.061 ± 0.005 


3.14 ±0.04 


31.3 ±2.4 


782.28 ± 0.14 


K5 


0.62 ±0.01 


0.838 ± 0.006 


0.27 ± 0.04 


0.298 ± 0.126 


3.13 ±0.07 


10.4 ±6.5 


782.41 ± 0.07 


K6 


0.58 ± 0.04 


0.843 ± 0.004 


0.30 ± 0.06 


0.061 ± 0.003 


3.14 ±0.01 


33.7 ±4.8 


782.48 ± 0.13 


NS 


1.43 ± 0.04 


0.953 ± 0.003 


0.00 ± 0.01 


0.270 ± 0.039 


3.11 ±0.14 


73.1 ± 1.6 


781.80 ±0.11 



9foK+K- = (3.76 ±0.04 fit to 



+ 0.15 +1.16 \ p Tr 

08 syst -0.48 mod; Vjt; V 



S/oTT+Tr- = (-1-43 ± 0.01 fit _o:o6 syst _o:60 mod) GeV 

p _ 9foK+K- _ g n _|_ n 1 - +0.2 +0..3 
■"-/o "~ — 2 ~ "'^ ± (J- J- fit -0.1 syst -3.9 mod 

94>foi = (2-78 toxil fit -0.05 syst + 1-31 mod) GeV ^ 



The first three quantities are the parameters directly ex- 
tracted from the fit while the other two are derived. The 
94, /o7 coupling is obtained using the formula 



9<j>Soi 



2M^ 



r^3BR((^-* 57^7707707). 



4.4.2 No Structure model 

In the No Structure model the resonant term is described 
with a single narrow meson pole, the /o(980), added to 
a continuum (f) — > w^Tr^'f background described by the 
three free parameters ao, ai and 61. The fit quality, hav- 
ing a P(x^) = 0.042, is a little worse than the best KL 
result but still acceptable (see Fig. 11). The different com- 
ponents of the fit are shown in Fig. 12. right. Again, for 
Mt^tt > 700 MeV the scalar contribution is clearly domi- 
nant. However, contrary to the KL case, the interference 



term gets negative in this region, so that the scalar term 
is slightly increased. The VMD fit parameters are listed 
in Tab. 5 while for the scalar term we have: 

Mf^ = (984.7 ± 0.4 fit t^-j syst) MeV 
9foK+K- = (0.40 ±0.04 fit 1°:^^ ,y,t) GeV 
9fo.+.- = (1.31 ±0.01 fit t^s syst) GeV 

Rfo = %^ = 0.09 ± 0.02 fit l^ots syst 

= (2.61 ± 0.02fit t"ofs syst) GeV-1 



ao= 4.19 ±0.01 fit tji^ syst 



ai = 1.31 ± 0.01 fit to.34 syst 



61 = -1.50 ±0.02fittn^ syst 

Here the only parameter that is not directly extracted 
from the fit is R/q. 

4.4.3 KL, NS combined results 

In Fig. 13 the comparison of the scalar contributions ob- 
tained from the six KL fit results shows stable result- 
ing shapes. On the contrary, comparing the KL and NS 
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Fig. 11. Fit result for No Structure model. Black dots are data while the solid line represents the NS resulting shape. 



curves there are differences at a level of few % (Fig. 13), 
mainly due to the interference term. In the same figure, 
the result of the fit obtained with the old /o(980)+cr(600) 
parametrization on 2000 data [8] is shown. The small bump 
below 500 MeV is now described by the 6*7- VMD inter- 
ference. By integrating the KL, NS distributions and nor- 
malizing to the (j) production cross section, an effective BR 
for the (j) Sj ^ 7r°7r°7 process is extracted: 

BR{(t> ^ 5*7 ^ Tr°Tr°-f) = 



^1 n7+0-01 +0.04 +0.05 \ „ 
yi.VI _o.03 fit _o.02 syst -O.O6 modj ^ 



The central value is given by the KL model with the best 
P(x^), the fit error has been evaluated as the maximum 
excursion obtained when varying the fit errors by ±la and 
the model error corresponds to the maximum variation of 
the central value with respect to the other five accepted 
fit results of the KL model and to the NS description. 
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Fig. 12. Resulting 7r°7r°7 contributions for Kaon Loop (left) and No Structure (right) models. 

4.4.4 Extrapolation to other ^/s points 



As a last check, we extrapolate both KL and NS fit results 
to the four closest ^/s points: 1019.55 MeV (42 ph'^), 
1019.65 MeV (77 pb"!), 1019.85 MeV (100 pb'^) and 
1019.95 MeV (15 pb~^), scaling for the integrated lumi- 
nosity. As shown in Figs. 14 and 15, a good agreement is 
obtained in all the M„„ and distributions. 




200 300 400 500 600 700 800 900 1000 

(MeV) 



Kaon Loop 
No Structure 
KLOE 2000 data 



4.5 Comparison with tt+tt 7 final state 



200 300 400 500 600 700 800 900 1000 

(MeV) 



In principle it is possible to compare these results with 

Fig. 13. M., distribution of the ^ 57 ^ 7r°7r°7 process ^^^^^ obtained by KLOE on the charged channel e+e" ^ 

obtained from the fit. Top: comparison between the six ac- + _ rooi xr • j. u i 

^ ^ TT^TT 7 [32J. However m this case the non resonant back- 

cepted results of the KL model. Bottom: comparison between , , , . , , , ,. ,. , .in . , 

ground (dommated by the p radiative tail) is much more 

the best KL fit, the NS resulting shape and the previous KLOE 

important, so that the extraction of the /o(980) signal 

measurement [8]. 

is more difficult and there is a reduced sensitivity to the 
presence of the cr(600). 
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Fig. 14. Dalitz plot projections for runs with different ^/s: 1019.55 MeV (top-left), 1019.65 MeV (top-right), 1019.85 MeV 
(bottom-left) and 1019.95 MeV (bottom-right). Data are reported in dots while the solid line represents the expected shape 
extrapolated from fit Kl of the Kaon Loop model. 



The results of the fit with the KL or NS model yield 
consistent values for the /o(980) mass and for the branch- 
ing ratio while large discrepancies are observed on the dif- 
ferential cross section, especially for the NS model. 

More precisely, by scaling by a factor two the extracted 
cross-section for the 7r°7r°7 case, we estimate BR(0 — ^ 
5*7 -> 7r+7r~7) = (2.14 ± 0.14) x lO^"*, to be compared 



to the values from 7r+7r-7 data: (2.1 ± 0.4) x lO""* (KL 
model) and (2.4 ± 0.5) x 10^"* (NS model). 

For the KL model, a still acceptable agreement on 
the /o(980) coupling constants is observed, although the 
(t(600) is not needed in the 7r^"7r^7 case. This is mainly 
due to the improved parametrization of the KL which had 
not been used for the charged case. For the NS model in- 
stead large differences are found on all couplings. 
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Fig. 15. Dalitz plot projections for runs witii different ^/s: 
(bottom-left) and 1019.95 MeV (bottom-right). Data are repi 
extrapolated from NS fit. 

5 Conclusions 

The analysis of the 7r'^7r'^7 final state presented here treats 
equally the two main production mechanisms, the VMD 
and 57 processes. The high statistics (145 pb^^ in a sin- 
gle y/s bin) allowed us to fit with two different theoretical 
models the Dalitz plot distribution. In the Kaon Loop 
model, the two low mass scalars /o(980) and a (600) are 



19.55 MeV (top-left), 1019.65 MeV (top-right), 1019.85 MeV 
;ed in dots while the solid line represents the expected shape 

required to adequately fit the data. In the case of the No 
Structure model, we instead find an acceptable fit with 
the /o(980) meson alone. In the latter formulation the low 
mass Mtt^^- behaviour is described by the three free parame- 
ters ao, ai and bi representing the continuum background. 
For both models, the resulting fit curve reproduces also 
the mass spectrum of all other ^/s bins around Af^. 
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ever, the physical interpretation is more difficult due to 
the presence of the continuum background which differs 
substantially in the tt^tt" and tt+tt" cases. 



A DifTerential cross section at y/s ~ 



The double differential 7r°7r°7 cross section can be written 
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where D^{s) is the ^ inverse propagator and g^^ is the 
coupling of the (j) to e+e". The general expression for a 



vector meson V is gv-y = y/3M^rvB{V e+e~)/a. The 
VMD parametrization contains two terms due to the ex- 
change of identical pions (M^^ vs M^^) and their inter- 
ference term. The full expression of the three coefficients 

Fi{m,m'^^) is the following: 

Fi(m2,m2 ) ^ ml^+2m^mtml-4mlml^ (7) 



+2 m^m^^ - 4 m^^m^^ + - 2 m^^s 



o22 2 ,r>4 2 r,24 

-2 m m^^m^^s + 2 m^^m^^s -2m m^^s 



F2(m2, ml^) = ml^ - m^m% -r^m^ml^ml^ (8) 



+6 m%^m%^ 2 m^m^^ - 4 m\^m%^ + m^^ 



+m m^i^ .s - 2 s + 2 m m^^s — 4m m^^m^^s 
+2 m^o"^^7* - m^m^^s + 2 m^^m^^s - 2 m^^s 

42 222, o2 22, 42 

mm ,m^^) — — ^ (9) 



The quantity Gp^u;{s,M^ ) is given by: 



, X r 



(10) 



The first term in the parenthesis is the only resonant com- 
ponent and includes the (jxv interference phase, which is 
set to 163° [36,37], and all the couplings involved in the 
reaction. Cp.^/^^ are complex coefficients that include the 



uncertainty arising from the recurrences of the p and w 
mesons and Sb^ is the phase of the amplitude when the p 
is the intermediate state. 
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